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The bright star e Aur (7 Aur, HD 31964, HR 1605; V max = 3 m 0; F01a+?) is an unusual 
eclipsing binary with a very long orbital period of 27.1 years (see Guinan & Dewarf 2002 
for a recent review). Its primary eclipse started in the summer 2009 and has naturally 
attracted the interest of many astronomers all over the world. The aim of this paper is 
to present our analysis of an extensive collection of archival and new photometry, and 
radial velocities (RVs), and provide a new, more precise ephemeris and orbital solution 
for the prediction of the current and future primary eclipses and a (not yet observed) 
secondary eclipse. Just prior to submission of this paper, Stefanik et al. (2010; hereafter 
ST) published their analysis of a comparable dataset for this same star. ST presented a 
new orbital solution and improved ephemeris for the binary but because the data anal- 
ysis approach presented here is significantly different and may provide a more accurate 
ephemeris, we have proceeded to publish our results also. 

We compiled and digitized a large collection of RVs from the literature, including ST's 
dataset of 515 RVs obtained at the Harvard-Smithsonian Center for Astrophysics (CfA). 
These data were augmented by our new series of electronic spectra from the Dominion 
Astrophysical Observatory (DAO) and the Ondfejov Observatory. Altogether, these RVs 
span an interval of 110 years. These RV observations are summarized in Table 1 and 
are plotted vs. time in Figure l 1 . We also collected and digitized light curves from all 
six previously observed eclipses. Additionally, for the 2010 eclipse, we used standard 
photoelectric V photometry obtained by PC, HB, DR, DS and MW at Hvar Observatory, 
CCD V-band photometry obtained by ML at the Hradec Kralove Observatory, and visual 
observations by PD reduced to Johnson V magnitude. These observations are listed in 
Table 2 and the individual eclipses are plotted in Figure 2 2 . 

1 RVs obtained during primary eclipse were not used in our solution because they are known to deviate from purely 
orbital motion. These eclipse RVs are not included in Table 1 or Figure 1. 

2 Some observations were not included because of their large scatter and/or unsufHcient coverage of a particular eclipse. 
We have also omitted extended datasets outside eclipse. These omitted data do not appear in either Table 2 or Figure 2. 



IBVS 5937 



Table 1. Journal of available RVs. 



from years 


observatory 


No. 


reference 


1899-1932 


Yerkes 


298 


Frost et al. (1929) 


1901-1913 


Post dam 


173 


Ludendorff (1924) 


1928-1958 


Mt.Wilson 


53 


Struve et al. (1958)* 


1970-1971 


Haute Provence 


18 


CasteUi (1977) 


1989-2009 


CfA 


515 


Stefanik et al. (2010) 


1994-2009 


DAO 


99 


this paper** 


2006-2009 


Ondfejov 


109 


this paper** 



RVs computed from the mean of 6 lines — Fe II 4123, Mg II 4481, Fe II 4508, Fc II 4515, Fe II 4576 and Fe II 4629 A. 
** RVs computed from the mean of 5 lines — Si II 6347, Si II 6371, Fc II 6417, Fe II 6433 and Fc II 6456 A. 

Table 2. Journal of photometric observations during primary eclipses. 



mid-eclipse 


observer 


passband* 


No. 


reference 


1848 


J.F.J.Schmidt 


pv 


39 


Ludendorff (1912) 


1875 


J.F.J.Schmidt 


pv 


69 


Ludendorff (1912) 


1902 


J.Plassmann 


pv 


29 


Ludendorff (1903) 


1902 


F.Schwab 


pv 


38 


Ludendorff (1903) 


1929 


C.M.Huffer & J.Stebbins 


pc 


98 


Huffer (1932) 


1956 


K.Gyldenkerne 


V 


131 


Gyldenkerne (1970) 


1956 


G . L ar sson- Leander 


V 


106 


Larsson-Leander (1959) 


1983 


J.L.Hopkins 


V 


130 


Schmidtke (1985) 


1983 


S.Ingvarsson 


V 


119 


Schmidtke (1985) 


2010 


Hvar Obs. 


V 


100 


this paper 


2010 


M.Lehky 


V 


21 


this paper 


2010 


P.Dubovsky 


pv(V) 


28 


this paper 



* Abbreviations 'pv' and 'pe' stand for photovisual and photoelectric, respectively. 



Here we provide more details of the new datasets. The DAO CCD spectra were ob- 
tained by SY and PDB and have a linear dispersion of 10 A mm -1 . The Ondfejov CCD 
spectra were obtained by PH, PS, MS, MW and a few other observers and have a dis- 
persion of 17 A mm -1 . Both the DAO and Ondfejov datasets cover the spectral region 
around 6300-6700 A. Their initial reductions were carried by SY and MS in IRAF. Rec- 
tification and RV measurements of the spectra were carried out by PC using the SPEFO 
(Horn et al. 1996, Skoda 1996) program's capability to compare direct and inverted line 
profiles. The zero point of the RV scale was determined by measurement of selected tel- 
luric lines (Horn et al. 1996). The Hvar dataset is actually UBV photometry carefully 
reduced to the standard system (Harmanec, Horn and Juza 1994). The Hradec Kralove 
CCD BVRI photometry was obtained with a 2.8/29 Pentacon auto lens and SBIG ST-5C 
CCD camera. The visual estimates by PD, reduced to Johnson y-band magnitude scale, 
were carried out using a modified version of Argelander's method developed by S. Otero 
(Stefl et al. 2003). It is based on a cone vision and calibration technique used to min- 
imize the effects of extinction and colour differences. We are making all new RVs and 
photometric datasets available with the electronic version of this paper 3 ; the remaining 
RV and photometric data are already accessible from the electronic version of ST. 



3 5937-tl.txt - t5.txt 
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Figure 1. RVs used in the orbital solution and the derived PHOEBE fit (curve). The vertical lines 
denote individual eclipses (during which RVs were not used in the solution). Plotted RVs are corrected 

for their individual 7 velocities. 



We used two independent programs, PHOEBE 0.31 (Prsa & Zwitter 2005) and FOTEL 
(Hadrava 2004), to derive new orbital solutions and formal light-curve solutions. All data 
sets were assigned weights inversely proportional to the squares of their rms errors derived 
from preliminary solutions. In FOTEL, we allowed calculations of individual systemic (7) 
velocities for individual spectrographs. Since PHOEBE can treat only a single RV set, 
we used RVs with individual 7 velocities subtracted. It turned out that the rms errors 
per observation for the RV sets in Table 1 were between 4 and 6 kms -1 . This indicates 
that the scatter is dominated by the intrinsic variations of the F star because the actual 
measurement errors are typically less than 1 kms" 1 . The RV solutions were used to derive 
the orbital eccentricity (e), longitude of periastron (u) and RV semiamplitude of a primary 
Ki, and the resulting values were then held fixed in the light-curve (LC) solutions. This 
is because the photometric data used only covers orbital phases near primary eclipse and, 
therefore, these data do not constrain the eccentric orbit. LC solutions were used to derive 
an improved ephemeris, assuming a mass ratio fixed at unity, and inclination fixed at 87°. 
The derived photometric period was held fixed again for the final iteration of the orbital 
solution, evaluated using the unconstrained system option in PHOEBE. 

The final photometric ephemerides (based exclusively on the LC solutions) are: 
Tprim.min. = HJD (2455402.8 ± 1.0) + (9890 d 26 ± d 62) x E (PHOEBE) , 
T prim . mill . = HJD (2455403.7 ± 1.1) + (9890 d 98 ± d 50) x E (FOTEL) . 

The epoch of primary minimum was allowed to vary independently for both the RV 
and LC solutions. We strongly prefer the more accurate value from photometry. For 
instance, the epoch of the primary minimum from the final RV solution in FOTEL at 
HJD 2455347 differs significantly from the above ephemerides. ST arrived at the same 
conclusions from their orbital solution; they obtained the epoch of the primary minimum 
at JD 2455136 (compared to their photometric minimum at JD 2455413). ST suggested 
that the gravitating companion responsible for the orbital motion need not be the same 
as the extended gaseous structure responsible for the eclipses. However, they also noted 
that intrinsic radial velocity variations in the F supergiant's atmosphere might bias the 
orbital solution, thereby accounting for the discrepancy between the photometric and RV 
solutions. We carried out an orbital solution in which the epoch of photometric mid- 
eclipse was held fixed and found that the resulting rms error was virtually identical to 
that of a solution converged with the epoch free to vary. This result strongly suggests 
that the discrepancy is due to intrinsic RV variations of the F supergiant and not due to 
asymmetry in the companion's structure. 
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Figure 2. Light curves from the last 6 eclipses, the current 2010 eclipse, and the PHOEBE fit (solid 

curve) are shown. Each measurement set is corrected to its individual 'zero level' magnitude. 
Mid-eclipse epochs have been centered using the new ephemeris and have been plotted on the same 

magnitude scale to facilitate visual comparison. 
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Table 3. New RV and LC solutions compared to those of Wright and Stefanik. 



element 


PHOEBE 


FOTEL 


Wright 


Stefanik 


T 

- 1 pcnastron 


2454596 ± 23* 


2454622 ± 97* 


2453130 ± 280* # 


2454515 ± 80* # 


T ■ 

1 pnm.mm. 


2455402.8 ± 1.0 


2455403.7 ± 1.1 


2455323 # 


2455413.8 ± 4.8 


T 

- L sec. mm. 


2451681 ± 120*+ 


2451610 ± 180*+ 


- 


- 


P(d) 


9890.26 ± 0.62 


9890.98 ± 0.50 


9890 (assumed) 


9896.0 ± 1.6 


e 


0.256 ± 0.012 


0.249 ± 0.015 


0.200 ± 0.034 


0.227 ± 0.011 


u; n 


41.2 ± 3.1 


43.3 ± 4.0 


346±11 


39.2 ± 3.4 


K x (kms- 1 ) 


14.40 ± 0.38 


14.30 ± 0.25 


15.00 ± 0.58 


13.84 ± 0.23 



* Errors from RV solutions. + Errors are semianalytical estimates. # Epochs recalculated for the authors' original periods. 
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Figure 3. A phase plot of all RVs used in the or- 
bital solution, and the derived PHOEBE fit (solid 
curve). Plotted RVs have been corrected for their 
individual 7 velocities. Note that PHOEBE also 
accounts for rotational effect during eclipse. 
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Figure 4. A phase plot of all photometric obser- 
vations used in the ephemeris calculation. Each 
measurement set has been corrected for individ- 
ual 'zero level' magnitude. 



We present our orbital solutions in Table 3, along with those of Wright (1970) and 
ST. A phase plot, using our new ephemeris from PHOEBE, of all RVs and photometry is 
shown in Figure 3. Note that our new solutions, obtained with two independent programs, 
agree within their respective errors. Our ephemerides predict the next primary mid-eclipse 
will occur on July 25-26, 2010, and the next secondary mid-eclipse in 2027. The previous 
secondary eclipse should have occurred in 2000. When compared to Wright, we obtain 
a significantly different orientation of the orbit in space (longitude of periastron u), a 
higher eccentricity (e), and a different epoch of the primary minimum. Our results are 
much closer to the ST solution, but we still disagree with ST by more than the estimated 
errors. At the request of the referee, we mention that the resulting relative photometric 
radii from the LC solutions were 0.045 and 0.216 from PHOEBE, and 0.058 and 0.218 
from FOTEL. We caution the reader, however, not to give these values much weight since 
neither program can treat disks; both assume two stellar bodies. 

Two important conclusions about e Aur, which disagree with the generally accepted 
model, follow from our study: 

1. Inspection of Figs. 2 and 3 shows that the idea of a central brightening inside 
the eclipse, interpreted as evidence of a hole in the disk (see, e.g., Carroll et al. 1991), 
should be reconsidered. Note that the 'flat' part of each recorded eclipse is different and 
what is seen are most probably the physical light variations, similar to the out-of-eclipse 
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variability. Of course, the final conclusion will come from a detailed analysis of colour 
changes and other types of observations and from the photometry secured this summer. 

2. The right panel of Figure 3 shows that claims of variability in the width and duration 
of individual observed eclipses, which have been used to infer a decline in the primary's 
radius over time (Saito 1986), are not supported by the data. It is apparent that the 
cyclic but irregular physical light variations affected the different eclipses differently. It 
will be difficult to obtain a 'pure' eclipsing light curve without a better understanding 
and quantitative description of these light changes. 
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